Origin of half-semimetallicity induced at interfaces of C-BN heterostructures 
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First-principles density functional calculations are performed in C-BN heterojunctions. It is shown 
that the magnetism of the edge states in zigzag shaped graphene strips and polarity effects in BN 
strips team up to give a spin asymmetric screening that induces half-semimetallicity at the interface, 
with a gap of at least a few tenths of eV for one spin orientation and a tiny gap of hundredths of 
eV for the other. The dependence with ribbon widths is discussed, showing that a range of ribbon 
widths is required to obtain half-semimetallicity. These results open new routes for tuning electronic 
properties at nanointerfaces and exploring new physical effects similar to those observed at oxide 
interfaces, in lower dimensions. 



A number of new physical phenomena have been dis- 
covered in the last few years at the interfaces between 
very different materials. The electronic reconstruction 
induced at these boundaries can give rise to metallic 
statesjy, magnetism [H, or even superconductivity [3], al- 
though the parent compounds were originally insulat- 
ing oxides. But not only boundaries of bulk materials 
are important: in two dimensional graphene nanoribbons 
(GNR) edges have become relevant, with peculiar elec- 
tronic states localized at the boundary of the ribbons 

Since its first experimental realization^, graphene has 
emerged as a prominent candidate to replace silicon in 
the development of high performance electronic devices. 
However, the extremely high mobility of charge carriers 
in graphene (ten times higher than that in silicon wafers 
used in microprocessors) poses a difficulty for the fabri- 
cation of nanodevices based on this material. A possible 
route to harness the conducting charges in graphene tran- 
sistors is the fabrication of graphene nanoribbons, where 
the lateral confinement of charge induces the opening of 
a gap that depends on the orientation of the GNR, and 
is inversely proportional to its width@,l. Interestingly, 
ribbons with zigzag shaped border (ZGNR) show a local 
ferrimagnetic structure at the edgeQ, and upon appli- 
cation of an in-plane external electric field can be tuned 
into a half-metallic state, with a gap opening for one 
spin component, and a metallic behaviour for the other, 
hence giving a full spin polarization of the conducting 
electrons 0. Theoretical calculations have shown that 
similar magnetic edge states exist in zigzag edged BN 
nanoribbons (ZBNNR)[l(], 11L The recent experimental 



realization of BN nanosheets 7(1 13] , opens new prospects 
for the combination of these two isostructural materials. 
Hence, a proper characterization of these hypothetical 
semimetal-insulator junctions is much needed. 

Different C/BN heterostructures have been proposed 
in the past, and substantial efforts have been devoted to 
the growth of composite sheets and nanotubesfljj]. Laser 
vaporization grown C-BN single walled nanotubes show 
traces of the patterning of segregated BN nanodomains 



embedded in the carbon network sequentially along the 
tube axis [25} . Recently, Du et. al. used first principles 
Molecular Dynamics calculations to show that hybrid C- 
BN nanotubes can be spontaneously formed via the con- 
nection of BNNR and GNR at room temperature [3], 
and Ding et. al. investigated the stability of C doped 
BNNR, showing that, under suitable conditions, GNR 
can be grown embedded in BN sheets 12|. Furthermore, 
it is reported that half- met allicity originates for suffi- 
ciently wide ZGNRs 12]. Tuning half-metallicity by edge 
modification in ZGNR and ZBNNR has been shown by 
others [ill 13 -HI, and understood in terms of a poten- 
tial difference between the two edges through chemical 
modification with electron accepting or donating groups. 
If this was the explanation for the half-metallicity in 
GNR embedded in BNNR, then there should be half- 
metallicity independently of the BNNR width. This 
doesn't seem to be the the results in reference [14[ 

suggest that there is a critical thickness for the BNNR to 
induce half-metallicity in GNRs. 

Here, I argue that the electronic properties of C+BN 
heterostructures can be described in terms of the ID 
equivalence of bidimensional oxide heterojunctions made 
of a metal (semimetal graphene) and an insulator (BN 
sheet). In particular, it will be shown that the polar- 
ity discontinuity in BN nanoribbon and the screening by 
the mobile GNR electrons are behind the electronic re- 
construction at the interface. Hence, a dependence with 
the BNNR width (polarity) and GNR width (number of 
available screening carriers) is expected. 

Ab initio pseudopotential density functional calcu- 
lations are performed in zigzag edged Graphene-BN 
superlattices, within the spin-polarized generalized- 
gradient approximation as implemented in the SIESTA 
method 0. Troullier- Martin type pseudopotentials fifij] 
and numerical atomic orbitals with double-^ plus polar- 
ization are used to describe the electronic valence states. 
The atomic positions are determined with a structural 
relaxation until the forces are smaller than 0.02eV/. A to- 
tal of 67 k-points are used to sample the Brillouin zone. 
Following the conventional nomenclature n-ZGNR and 
771-ZBNNR combine to give a (n, m)-superlattice with n 
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FIG. 1: (a) Diagram of a (5,9) C-BN superlattice with 5 and 
9 zigzag chains for the C-ribbon and BN-ribbon, respectively. 
White, black, and dark gray circles represent C, N and B 
atoms. The dashed line box shows the cell of repetition un- 
der periodic boundary conditions. Calculated C-C, B-N, C- 
B and C-N bond lengths are also shown, (b) Profile of the 
macroscopic electronic effective potential averaged over the 
graphene plane. The origin of the potential was shifted to the 
N-edge of the superlattice. 



zigzag chains of graphene and m zigzag chains of BN (Fig. 
la). A set of systems with 2 < n < 12 and 2 < m < 11 
are considered, with superlattice's periodicities in the 
range 1.7-4.7nm. The smaller systems (n + m = 8) show 
non-magnetic (NM) semiconducting properties, with a 
direct energy gap that decreases with the width of the 
graphene ribbon. For wider systems however the behav- 
ior becomes semimetallic, with edge states that give rise 
to a narrow band near the Fermi level, enabling the pos- 
sibility of magnetic orderings induced by the electron- 
electron interaction Indeed, spin-polarized edge states 
are obtained in the calculations, that correspond to an- 
tiferromagnetic (AF) alignments of the spin moments at 
opposite edges of the heterostructure (Fig. [21 and Sup- 
porting Information) , with the AF state being a few meV 
lower than NM, in agreement with calculations of mag- 
netic edge states in isolated ZGNR, and ZGNR embed- 
ded in ZBNNR 



12]. Double periodicity unit cells were 



used to explore other magnetic orderings but this AF- 
phase remained as the ground state. Notice that hybrid 
exchange-DFT functionals predict a further stabilization 
of the AF alignment in isolated ZGNR 21]. Notice also 
the half-semimetallic character of the AF state in Fig. 2, 
with an apparent gap for a spin, and two bands that seem 
to cross the Fermi level for /3 spin. The crossing is not 
allowed by symmetry because each sublattice of C atoms 
is linked either to B or N breaking the sub-lattice's equiv- 
alency and giving a tiny gap also for the /3-spin bands. 
Hence, the term half-semimetal is more appropriate than 
half-metal commonly used in the literature. 

The nature of these edge bands near the Brillouin zone 
boundary (X-point at k=7r) deserves further clarification. 
Zigzag GNRs have insulating magnetic edge states with 



parallel ordering at each edge, and antiparallel alignment 
between the two edges. On the other hand, zigzag BN- 
NRs have insulating, nonmagnetic edge states with unoc- 
cupied state localized at the B site and occupied valence 
state localized at the N side. These states can be dramat- 



ically affected by edge passivation [ll[. When GNR and 
BNNR are attached the mixing of tt orbitals of C, B and 
N at each edge gives rise to four sets of bands that cor- 
respond to the bonding and antibonding states between 
C-N and C-B (Fig. 3a). The relevant bands, close to 
the Fermi level, are the occupied bonding C-B (ttb), and 
unoccupied antibonding C-N (tt^), each of them being a 
localized state at carbon atoms close to the B or N edge, 
respectively. Figure 3.b-e presents evidences of the highly 
localized states at X by showing a color contour plots for 
the energy bands near the Fermi level projected on the C 
atomic orbitals at the atoms close to the N and B edges. 
The localization gradually decreases as we depart from 
the X-point. The bonding C-N and antibonding C-B lay 
~5eV below and above the Fermi level respectively. Ab 
initio calculations have shown that the bandgaps in both 
C and BN zigzag nanoribbons are inversely proportional 
to their width, and consequently the energy separation 
between ttb and in C-BN super lattices will depend 
on the width of each strip. Hence, a dependence of the 
half-semimetallic gap on the widths of the strips is antic- 
ipated, and will be discussed later. 

The different electronegativities of B and N make bare 
ZBNNR the two-dimensional equivalent to a polar slab, 
with a type 3 termination, where alternating charges are 
ordered perpendicularly to the edge, as opposed to arm- 
chair BNNR which have type 1 boundaries and alter- 
nating charges parallel to the interface 22| . Polarization 




FIG. 2: (color online) Electronic bandstructure of C-BN su- 
perlattices in the antiferromagnetic (AF) groundstate and 
non-magnetic (NM) state for the (5,9)-superlattice. Dashed 
(blue) and solid (red) lines correspond to a and ft spins. The 
presence of an energy gap, A, for the a spin is apparent in 
the AF state. 
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FIG. 3: (color online) (a) Energy diagram for localized states 
close to the C-N (left) and C-B (right) edges, for a (blue) 
and /3 (red) spins. GNR edge states have opposite spin char- 
acter at each side of the ribbon, and couple to the B and N 
orbitals. Near the Fermi level (dashed horizontal line), the rel- 
evant states are occupied ns and unoccupied 7r^, with some 
spin asymmetry coming from the original GNR spin edge po- 
larization. The weights on the relevant bands of the C atomic 
orbitals at the N and B edges are shown in b, c, d and e for 
each spin, with red/blue being a high/negligible contribution 
to the electronic wavefunction. The Fermi level is taken at 
and the bands are plotted along the T-X direction. 



lowers the potential felt by electrons at the nitrogen edge 
relative to the boron edge. Charge redistributions near 
the interface partially compensate this edge instability, 
particularly for wider ribbons, where the number of elec- 
trons that can participate in the screening is higher(23j|. 
It is then expected that mobile electrons coming from 
semimetallic graphene in contact with these ZBNNR will 
increase the screening. But this also means that there 
is a charge asymmetry in the GNR in contact with BN 
with the interface C-N being lower in energy than the 
C-B interface, and hence inducing a charge transfer to 
the N side. This charge transfer is confirmed by Mul- 
liken population analysis, with a slightly larger charge 
for C atoms close to the N interface, than for C atoms 
close to the B edge. There is also a minor charge in- 
crease for the B atom close to the B-edge 14], but its role 
in origin of half-metallicity is only partial. The situation 
very much resembles that of a GNR under the presence 
of a transverse electric field, except that the nature of 
the edge states is different, with occupied bands near the 
Fermi level localized at the same edge (B-side) for both 
spin polarizations. Atomic relaxations also contribute to 
the screening, with bonds becoming shorter at the N-side 
and longer at the B-side. 

Figure lb shows the macroscopic average [24] of the 
electronic effective potential in the plane of the (5,9)- 
superlattice. This potential profile shows a drop of 
~2.3eV between B and N edges, inducing an effective 
electric field of ~0.27V/A, which is at the brink of the 
critical electric field needed to induce half-metallicity 
in the isolated ZGNR|9]. The potential drop increases 
slightly with the BN widths studied here but the decrease 
of the effective polarization charge density ensures that it 
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FIG. 4: (color online) Direct band gaps for a- (blue up tri- 
angles) and /3-spins (red down triangles) as a function of the 
GNR width in the antiferromagnetic ground state for C+BN 
superlattices. The (green) arrow shows the dependence of the 
a-spin gap on the BNNR width (also displayed in the inset). 
Superlattices made from GNR thinner than 15 and BNNR 
thicker than 10 present a marked half-semimetallicity with a 
O.leV gap for a-spins and a tiny gap for /3-spins. 



will saturate for sheets wide enough. For wider graphene 
strips, the screening is asymmetric and mainly localized 
a few (4 or 5) zigzag chains from the B-edge, so that 
the effective electric field inside the graphene ribbon will 
be substantially screened far from this edge. This means 
that the relative widths of the GNR and BNNR strips can 
be used to tune the electronic and magnetic properties of 
these systems. The most important result of this work is 
shown in Fig. 4 where the band gaps for a and /3 spins are 
plotted as a function of the strip widths. In the limit of 
broad GNR, the semimetallic behavior is recovered, with 
zero gap for both spins. If the GNR is too thin, the AF 
state is unstable and the systems becomes NM and insu- 
lating. However, lattices with sufficiently wide BN strips 
(to > 6) and sufficiently narrow C ribbons (n < 8) can 
show half-semimetallic properties with a gap for a-spin 
states significantly larger (above O.leV) than for /3-spin 
(below 0.03eV). It is important to notice that the origin 
of this effect is not simply the chemical modification of 
the ZGNR by B or N doping, as the spin polarization 
disappears for narrow BN strips (inset in Fig. 4). 

As mentioned before, hybrid C-BN nanotubes have al- 
ready been synthesized experimentally [25}, and predic- 
tions have been made for the spontaneous formation of 
single-walled armchair nanotubes from the hybrid con- 
nection of BNNR and GNR 26] . Unzipping these hybrid 
nanotubes would be a route to fabricate C-BN superlat- 
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tices as the ones reported in this work. However, it can 
be shown [? ] that even in the tubular geometry the elec- 
tronic and magnetic properties described above can be 
obtained, as long as there is a zigzag edge and the rela- 
tive width of BN and C regions are conveniently tuned, 
following the gap width-dependence sketched in Fig. 4. 

In addition to the demonstration that unconventional 
intcrfacial effects are not limited only to planar junc- 
tions in epitaxial complex oxides heterostructures but 
are present in lower dimensions, synthesis of these C-BN 
nanostructure offers an interesting avenue for the design 
of carbon-based nanospintronic devices. A zigzag-edged 
nanostriction of C-BN in a graphene ribbon, for example, 
would open a gap and induce a high spin-polarization for 
the transmitted electrons without the need for applied 
external electric fields, paving the way to efficient spin- 
injection in carbon structures. Surely different combi- 
nation of materials can be conceived to fabricate bidi- 
mensional heterostructures and take advantage of these 
unusual physical phenomena that result from the com- 
bination of spin-polarized edge states and polarization 
induced charge dipoles at the edges. Further work on 
growing bidimensional superlattices that combine the ex- 
iting possibilities of nanosheets is encouraged. 
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First-principles density functional calculations are performed in C-BN heterojunctions. It is shown 
that the magnetism of the edge states in zigzag shaped graphene strips and polarity effects in BN 
strips team up to give a spin asymmetric screening that induces half-semimetallicity at the interface, 
with a gap of at least a few tenths of eV for one spin orientation and a tiny gap of hundredths of 
eV for the other. The dependence with ribbon widths is discussed, showing that a range of ribbon 
widths is required to obtain half-semimetallicity. These results open new routes for tuning electronic 
properties at nanointerfaces and exploring new physical effects similar to those observed at oxide 
interfaces, in lower dimensions. 



A number of new physical phenomena have been dis- 
covered in the last few years at the interfaces between 
very different materials. The electronic reconstruction 
induced at these boundaries can give rise to metallic 
statesjy, magnetism or even superconductivity [3], al- 
though the parent compounds were originally insulat- 
ing oxides. But not only boundaries of bulk materials 
are important: in two dimensional graphene nanoribbons 
(GNR) edges have become relevant, with peculiar elec- 
tronic states localized at the boundary of the ribbons 

Since its first experimental realization^, graphene has 
emerged as a prominent candidate to replace silicon in 
the development of high performance electronic devices. 
However, the extremely high mobility of charge carriers 
in graphene (ten times higher than that in silicon wafers 
used in microprocessors) poses a difficulty for the fabri- 
cation of nanodevices based on this material. A possible 
route to harness the conducting charges in graphene tran- 
sistors is the fabrication of graphene nanoribbons, where 
the lateral confinement of charge induces the opening of 
a gap that depends on the orientation of the GNR, and 
is inversely proportional to its width^, 0]. Interestingly, 
ribbons with zigzag shaped border (ZGNR) present a lo- 
cal ferrimagnetic structure at the edgeQ, and upon appli- 
cation of an in-plane external electric field can be tuned 
into a half-metallic state, with a gap opening for one 
spin component, and a metallic behavior for the other, 
hence giving a full spin polarization of the conducting 
electrons Theoretical calculations have shown that 
similar magnetic edge states exist in zigzag edged BN 
nanoribbons (ZBNNR)@, [l(L_ The recent experimental 
realization of BN nanosheets UJ, , opens new prospects 
for the combination of these two isostructural materials. 
Hence, a proper characterization of these hypothetical 
semimetal-insulator junctions is much needed. 

Different C/BN heterostructures have been proposed 
in the past, and substantial efforts have been devoted to 
the growth of composite sheets and nanotubesflii]. Laser 
vaporization grown C-BN single walled nanotubes show 
traces of the patterning of segregated BN nanodomains 



embedded in the carbon network sequentially along the 
tube axis [3]. Recently, Du et. al. used first principles 
Molecular Dynamics calculations to prove that hybrid C- 
BN nanotubes can be spontaneously formed via the con- 
nection of BNNR and GNR at room temperature [la], and 
Ding et. al. investigated the stability of C doped BNNR, 
showing that, under suitable conditions, GNR can be 
grown embedded in BN sheets [IB]. Furthermore, it is 
reported that half-metallicity originates for sufficiently 
wide ZGNRsflfl- Tuning half-metallicity by edge mod- 
ification in ZGNR and ZBNNR has been demonstrated 
by others 1(J 11 - 3] , and understood in terms of a po- 
tential difference between the two edges through chemical 
modification with electron accepting or donating groups. 
If this was the explanation for the half-metallicity in 
GNR embedded in BNNR, then there should be half- 
metallicity independently of the BNNR width. This 
doesn't seem to be the the results in reference [3 

suggest that there is a critical thickness for the BNNR to 
induce half-metallicity in GNRs. 

Here, I argue that the electronic properties of C+BN 
heterostructures can be described in terms of the ID 
equivalence of bidimensional oxide heterojunctions made 
of a metal (semimetal graphene) and an insulator (BN 
sheet). In particular, it will be shown that the polarity 
discontinuity in BN nanoribbon and the screening by the 
mobile electrons in graphene are behind the electronic re- 
construction at the interface. Hence, a dependence with 
the BNNR width (polarity) and GNR width (number of 
available screening carriers) is expected. 

Ab initio pseudopotential density functional calcu- 
lations are performed in zigzag edged Graphene-BN 
superlattices, within the spin-polarized generalized- 
gradient approximation as implemented in the SIESTA 
method [20(. Troullier- Martin type pseudopotentials 21 
and numerical atomic orbitals with double-^ plus polar- 
ization are used to describe the electronic valence states. 
The atomic positions are determined with a structural 
relaxation until the forces are smaller than 0.02eV/. A to- 
tal of 67 k-points are used to sample the Brillouin zone. 
Following the conventional nomenclature n-ZGNR and 
TO-ZBNNR combine to give a (n, m)-superlattice with n 
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FIG. 1: (a) Diagram of a (5,9) C-BN superlattice with 5 and 
9 zigzag chains for the C-ribbon and BN-ribbon, respectively. 
White, black, and dark gray circles represent C, N and B 
atoms. The dashed line box shows the cell of repetition un- 
der periodic boundary conditions. Calculated C-C, B-N, C- 
B and C-N bond lengths are also shown, (b) Profile of the 
macroscopic electronic effective potential averaged over the 
graphene plane. The origin of the potential was shifted to the 
N-edge of the superlattice. 



zigzag chains of graphene and m zigzag chains of BN 
(Fig. la). A set of systems with 2 < n < 12 and 
2 < m < 11 are considered, with superlattice's peri- 
odicities in the range 1.7-4.7nm. The smaller systems 
(n + m = 8) show non-magnetic (NM) semiconducting 
properties, with a direct energy gap that decreases with 
the increase in the width of the graphene ribbon. For 
wider systems however the gap is reduced and the sys- 
tem semimetallic, with edge states that give rise to a nar- 
row band near the Fermi level, enabling the possibility 
of magnetic orderings induced by the electron-electron 
interaction^. Indeed, spin-polarized edge states are ob- 
tained in the calculations, that correspond to antifcrro- 
magnetic (AF) alignments of the spin moments at op- 
posite edges of the heterostructure (Fig. [2]), with the 
AF state being a few meV lower than NM, in agree- 
ment with calculations of magnetic edge states in isolated 
ZGNR, and ZGNR embedded in ZBNNrQ. Double pe- 
riodicity unit cells were used to explore other magnetic 
orderings observed for isolated BNNR[9], but the AF- 
phase remained as the ground state. Notice that hybrid 
cxchangc-DFT functionals predict a further stabilization 
of the AF alignment in isolated ZGNR [H]. Notice also 
the half-semimctallic character of the AF state in Fig. 2, 
with an apparent gap for a spin, and two bands that seem 
to cross the Fermi level for /3 spin. The crossing is not 
allowed by symmetry because each sublattice of C atoms 
is linked either to B or N, breaking the sublattice's equiv- 
alency and giving a tiny gap also for the /3-spin bands. 
Hence, the term half-semimetal is more appropriate than 
half-metal commonly used in the literature. 

The nature of these edge bands near the Brillouin zone 
boundary (X-point at k=7r) deserves further clarifica- 



tion. Hydrogen passivated zigzag GNRs have insulat- 
ing magnetic edge states with parallel ordering at each 
edge, and antiparallel alignment between the two edges. 
On the other hand, zigzag BNNRs have insulating, non- 
magnetic edge states with unoccupied state localized at 
the B site and occupied valence state localized at the N 
side. These states can be dramatically affected by edge 
passivation [13, [13]. When GNR and BNNR are attached 
the mixing of it orbitals of C, B and N at each edge gives 
rise to four sets of bands that correspond to the bonding 
and antibonding states between C-N and C-B (Fig. 3a). 
The relevant bands, close to the Fermi level, are the oc- 
cupied bonding C-B (ttb), and unoccupied antibonding 
C-N (7r^), each of them being a localized state at car- 
bon atoms close to the B or N edge, respectively. Figure 
3.b-e presents evidences of the highly localized states at 
X by showing a color contour plots for the energy bands 
near the Fermi level projected on the C atomic orbitals 
at the atoms close to the N and B edges. The localization 
gradually decreases as we depart from the X-point. The 
bonding C-N and antibonding C-B lay ~5eV below and 
above the Fermi level respectively. Ab initio calculations 
have shown that the bandgaps in both C and BN zigzag 
nanoribbons are inversely proportional to their width, 
and consequently the energy separation between ttb and 
ir* N in C-BN superlattices will depend on the width of 
each strip. Hence, a dependence of the half-semimctallic 
gap on the widths of the strips is anticipated, and will be 
discussed later. 

The different electronegativities of B and N make bare 
ZBNNR the two-dimensional equivalent to a polar slab, 
with a type 3 termination, where alternating charges are 
ordered perpendicularly to the edge, as opposed to arm- 




FIG. 2: (color online) Electronic bandstructure of C-BN su- 
perlattices in the antiferromagnetic (AF) groundstate and 
non-magnetic (NM) state for the (5,9)-superlattice. Dashed 
(blue) and solid (red) lines correspond to a and /3 spins. The 
presence of an energy gap, A, for the a spin is apparent in 
the AF state. 
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FIG. 3: (color online) (a) Energy diagram for localized states 
close to the C-N (left) and C-B (right) edges, for a (blue) 
and /3 (red) spins. GNR edge states have opposite spin char- 
acter at each side of the ribbon, and couple to the B and N 
orbitals. Near the Fermi level (dashed horizontal line), the rel- 
evant states are occupied ns and unoccupied 7r^, with some 
spin asymmetry coming from the original GNR spin edge po- 
larization. The weights on the relevant bands of the C atomic 
orbitals at the N and B edges are shown in b, c, d and e for 
each spin, with red/blue being a high/negligible contribution 
to the electronic wavefunction. The Fermi level is taken at 
and the bands are plotted along the T-X direction. 



chair BNNR which have type 1 boundaries and alter- 
nating charges parallel to the interface [HI ■ Polarization 
lowers the potential felt by electrons at the nitrogen edge 
relative to the boron edge. Charge redistributions near 
the interface partially compensate this edge instability, 
particularly for wider ribbons, where the number of elec- 
trons that can participate in the screening is higher(25j|. 
Atomic relaxations also contribute to the screening, with 
bonds becoming shorter at the N-side and longer at the 
B-side. It is then expected that mobile electrons com- 
ing from semimetallic graphene in contact with these 
ZBNNR will increase the screening. But this also means 
that there is a charge asymmetry in the GNR in con- 
tact with BN with the interface C-N being lower in en- 
ergy than the C-B interface, and hence inducing a charge 
transfer to the N side. This charge transfer is confirmed 
by Mulliken population analysis, with a slightly larger 
charge for C atoms close to the N interface, than for C 
atoms close to the B edge. This result is at odds with the 
hypothesis that the origin of the half-mctallicity in the 
GNR is the charge transfer from the carbon to boron 
atoms [lg. Although the occupied states close to the 
Fermi level are mainly localized at the atoms close to 
this side (Fig. 3), the fact that for very narrow ZGNR 
half-semimetallicity is not observedjlq proves that this 
postulate is incorrect. 

The situation very much resembles that of a GNR un- 
der the presence of a transverse electric field, except that 
the nature of the edge states is different, with occupied 
bands near the Fermi level localized at the same edge 
(B-side) for both spin polarizations. Figure lb shows the 
macroscopic average [26] of the electronic effective poten- 
tial in the plane of the (5,9)-superlattice. This potential 



profile shows a drop of ~2.3eV between B and N edges, 
inducing an effective electric field of ~0.27V/A, which is 
at the brink of the critical electric field needed to induce 
half-metallicity in the isolated ZGNR0. The potential 
drop increases slightly with the BN widths studied here 
but the decrease of the effective polarization charge den- 
sity ensures that it will saturate for sheets wide enough. 
For wider graphene strips, the screening is asymmetric 
and mainly localized a few (4 or 5) zigzag chains from 
the B-edge, so that the effective electric field inside the 
graphene ribbon will be substantially screened far from 
this edge. This means that the widths of the GNR and 
BNNR strips can be used to tune the electronic and mag- 
netic properties of these systems. 

The most important result of this work is shown in 
Fig. 4 where the band gaps for a and /3 spins are plotted 
as a function of the strip widths. In the limit of broad 
GNR, the semimetallic behavior is recovered, with zero 
gap for both spins. If the GNR is too thin, the AF state 
is unstable and the systems becomes NM and insulat- 
ing. However, lattices with sufficiently wide BN strips 
(to > 6) and sufficiently narrow C ribbons (n < 8) can 
show half-semimetallic properties with a gap for a-spin 
states significantly larger (above O.leV) than for /3-spin 
(below 0.03eV). It is important to notice that the origin 
of this effect is not simply the chemical modification of 
the ZGNR by B or N doping, as the spin polarization 
disappears for narrow BN strips (inset in Fig. 4). 

As mentioned before, hybrid C-BN nanotubes have al- 
ready been synthesized experimentally Q, and predic- 
tions have been made for the spontaneous formation of 
single-walled armchair nanotubes from the hybrid con- 
nection of BNNR and GNR[l5j]. Unzipping these hybrid 
nanotubes would be a route to fabricate C-BN superlat- 
tices as the ones reported in this work. However, it can 
be shown that even in the tubular geometry the elec- 
tronic and magnetic properties described above can be 
obtained, as long as there is a zigzag edge and the widths 
of BN and C ribbons are conveniently tuned, following 
the gap width-dependence sketched in Fig. 4. Curva- 
ture does not play an essential role, and (n, n) zigzag 
nanotubes with n > 5 can show these half-metallic prop- 
erties. Details will be presented elsewhere 27 1. 

In addition to the demonstration that unconventional 
interfacial effects are not limited only to planar junc- 
tions in epitaxial complex oxides heterostructures but 
are present in lower dimensions, synthesis of these C-BN 
nanostructure offers an interesting avenue for the design 
of carbon-based nanospintronic devices. A zigzag-edged 
nanostriction of C-BN in a graphene ribbon, for example, 
would open a gap and induce a high spin-polarization for 
the transmitted electrons without the need for applied 
external electric fields, paving the way to efficient spin- 
injection in carbon structures. Surely different combi- 
nation of materials can be conceived to fabricate bidi- 
mensional heterostructures and take advantage of these 
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FIG. 4: (color online) Direct band gaps for ot- (blue up tri- 
angles) and /3-spins (red down triangles) as a function of the 
GNR width in the antiferromagnetic ground state for C4-BN 
superlattices. The (green) arrow shows the dependence of the 
a-spin gap on the BNNR width (also displayed in the inset). 
Superlattices made from GNR thinner than 15 and BNNR 
thicker than 10 present a marked half-semimetallicity with a 
O.leV gap for a-spins and a tiny gap for /3-spins. 



unusual physical phenomena that result from the com- 
bination of spin-polarized edge states and polarization 
induced charge dipoles at the edges. Further work on 
growing bidimensional superlattices that combine the ex- 
iting possibilities of nanosheets is encouraged. 

This work was supported by the CSIC through a JAE- 
Doc fellowship and PTE. 2009601025. Computational 
resources from the Theory and Simulation group are ac- 
knowledged. The author is thankful to P. Ordejon, E. 
Artacho and E. Canadell for fruitful discussions and com- 
ments. 



* miguel.pruneda@cin2.es 
[1] A. Ohtomo and H. Y. Hwang. Nature 427. 423 (2004). 
[2] A. Brinkman, M. Huijben, M. van Zalk, J. Huijben, 



U. Zeitler, J. C. Maan, W. G. van der Wiel, G. Rinjders, 
D. H. A. Blank, and J. Hilgenkamp, Nat. Materials 6, 
493 (2007). 

[3] N. Reyren, S. Thiel, A. D. Caviglia, L. F. Kourkoutis, 

G. Hammerl, C. Richter, C. W. Schneider, T. Kopp, A.- 

S. Ruetschi, D. Jaccard, et al., Science 317, 1196 (2007). 
[4] M. Fujita, K. Wakabayashi, K. Nakada, and K. J. Kusak- 

abe, J. Phys. Soc. Jpn. 65, 1920 (1996). 
[5] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
[6] M. Y. Han, B.Ozyilmaz, Y. Zhang, and P. Kim, Phys. 

Rev. Lett. 98, 206805 (2007). 
[7] Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 

Lett. 97, 216803 (2006). 
[8] Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 

347 (2006). 

[9] V. Barone and J. E. Peralta, Nano Letters 8, 2210 (2008). 
[10] F. Zheng, G. Zhou, Z. Liu, J. Wu, W. Duan, B.-L. Gu, 

and S. B. Zhang, Phys. Rev. B 78, 205415 (2008). 
[11] D. Pacile, J. C. Meyer, JC. O. Girit, and A. Zettl, Appl. 

Phys. Lett. 92, 133107 (2008). 
[12] C. Jin, F. Lin, K. Suenaga, and S. Iijima, Phys. Rev. 

Lett. 102, 195505 (2009). 
[13] M. Kawaguchi, Adv. Mater. 9, 615 (1997). 
[14] S. Enouz, O. Stephan, J.-L. Colliex, and A. Loiseau, 

Nano Letters 7, 1856 (2007). 
[15] A. Du, Y. Chen, Z. Zhu, G. Lu, and S. C. Smith, J. Am. 

Chem. Soc. 131, 1682 (2009). 
[16] Y. Ding, Y. Wang, and J. Ni, Appl. Phys. Lett. 95, 

123105 (2009). 

[17] E.-J. Kan, Z. Li, J. Yang, and J. G. Hou, J. Am. Chem. 
Soc. 130, 4224 (2008). 

[18] S. Dutta, A. K. Manna, and S. K. Pati, Phys. Rev. Lett. 
102, 096601 (2009). 

[19] M. Wu, X. Wu, Y. Gao, and X. C. Zeng, Appl. Phys. 
Lett. 94, 223111 (2009). 

[20] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal, J. Phys.: Con- 
dens. Matter 14, 2745 (2002). 

[21] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

[22] L. Pisani, J. A. Chan, B. Montanari, and N. M. Harrison, 

Phys. Rev. B 75, 064418 (2007). 
[23] G. Lee and K. Cho, Phys. Rev. B 79, 165440 (2009). 
[24] P. W. Tasker, J. Phys. C: Solid State Phys. 12, 4977 

(1979). 

[25] C.-H. Park and S. G. Louie, Nano Letters 8, 2200 (2008). 
[26] A. Baldereschi, S. Baroni, and R. Resta, Phys. Rev. Lett. 

61, 734 (1988). 
[27] J. M. Pruneda, to be published. 



